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ANALYSIS  OF  LAMINATED  COMPOSITE  PLATES  USING  HYBRID-STRESS  ELEMENTS 
Deng  Liangbo 

Summary 

Based  on  the  modified  complementary  energy  principle, 
the  paper  describes  a  hybrid-stress  rectangular  flexural 
element  suitable  for  characteristics  of  laminated  fiber 
composite  plates.  This  element  can  take  into  account  the 
effects  of  transverse  shear  deformation  and  local  warping. 

In  the  dynamic  analysis,  the  interior  element  displacements 
are  assumed  to  be  compatible  with  the  boundary  displace¬ 
ments.  The  element  has  the  characteristics  of  fewer 
degrees  of  freedom,  convenience  in  application,  and  high 
accuracy  in  numerical  calculations.  •  •••-"'■  •  " 

I.  Introduction 

When  applying  the  method  of  finite  elements  to  analyze  lamina¬ 
ted  plates  made  of  composites,  not  only  should  transverse  shear 
deformations  be  considered,  but  also  the  local  warping  effect.  In 
this  respect,  the  hybrid-stress  finite  elements  method  has  unique 
advantages.  When  constituting  elements  in  [1],  the  interlaminar 
displacement  of  nodal  points  was  taken  as  a  degree  of  freedom. 
Hence,  the  number  of  degrees  of  freedom  of  an  element  will  rise 
with  an  increase  in  the  number  of  layers,  thus  incurring  diffi¬ 
culties  in  numerical  calculations,  especially  in  multilaminated 
plates.  An  element  is  proposed  in  [2];  only  a  set  of  stress 
parameters  is  assumed  in  the  element  interior,  and  a  hypothesis 
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of  YNS  theory  (hypothesis  of  rectilinear  normal  line)  is  applied 
on  the  element  boundary,  therefore  the  local  warping  effect  can  not 
be  taken  into  account.  With  the  foregoing  in  mind,  the  author 
conceptualizes  the  title  element;  the  hypothesis  of  YNS  theory  is 
adopted  for  the  element  boundary,  and  stress  parameters  are  assumed 
layer  by  layer  in  the  element  interior.  Hence,  the  element  can  take 
into  account  both  the  effect  of  transverse  shearing  deformation  and 
the  local  warping  effect  with  fewer  degrees  of  freedom. 

II.  Theory 

In  the  analysis  of  laminated  plates,  the  modified  arbitrary 
function  of  the  rest  energy  is: 

y'ts'nio  *  >dv 

+  fsaai*<p’  (i)  }ds)J 

In  the  equation, 

{ c ( 1 5  }  is  the  stress  vector  of  the  i-th  subzone  (i-th  layer); 

{s (i) }  is  the  elasticity  matrix  of  the  i-th  subzone; 

{u (i) }  is  the  boundary  displacement  vector  of  the  i-th  subzone; 
{p^1^  }  is  the  boundary  force  vector  of  the  i-th  subzone; 

{p^  is  the  assigned  boundary  force  vector  of  the  i-th 
sub zone ; 

9vn^  is  all  boundary  3vni  =  Suni  +  Scrni  +  Sni  °f  the  i-th  subzone 
of  the  n-th  element; 

v  .  is  the  volume  of  the  i-th  subzone  of  the  n-th  element; 
ni 

Sa^  is  the  delimited  boundary  of  the  i-th  subzone  boundary 
force  of  the  n-th  element; 

Sun^  is  the  delimited  boundary  in  terms  of  displacements; 
and  is  the  boundary  adjacent  to  other  elements. 

When  one  applies  the  arbitrary  function  (1),  the  assumed 
stress  vector  should  satisfy  the  homogeneous  equilibrium  equation 


-I 


Uv, 


{P 


[u  }  ds  -r 


(1) 
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in  each  sabzone.  At  the  boundary  Su  .  ,  {u  ^  }  =  {u  ^  }  should  be 

/  •  \  m 

satisfied;  in  the  equation,  {u  1  }  is  the  known  displacement  vector 
at  the  boundary. 


We  assume  that  the  stress  field  of  each  layer  is: 

{cr (i>  >  =  [ A  ]{{3  li>  >  (2) 

In  the  equation,  [A ^ ]  is  the  coefficient  matrix  of  stress  param¬ 
eters;  this  is  a  function  of  x,  y  and  z.  }  is  the  stress 

parameter  of  the  i-th  layer. 


In  terms  of  the  stress  parameter,  the  unit  boundary  force  vec¬ 
tor  can  be  expressed  as: 

}  (3) 

In  the  equation,  [R  ^  ]  is  constituted  with  the  boundary  value  of 

the  matrix  [A^]  in  satisfying  the  boundary  force  equation  P^= 

=  a  .  .  n  .  . 

1 D  ] 

We  assume  that  the  boundary  displacement  of  the  i-th  layer  is: 

{u'i>  }  =  [L(i>  >  (4) 

In  the  equation,  {6^  is  the  nodal  point  displacement  (in  the 
broad  sense)  of  the  i-th  layer.  The  boundary  displacement  field 
(4)  should  satisfy  the  displacement  compatibility  conditions 
among  elements. 

If  we  substitute  Eqs.  (2),  (3)  and  (4)  into  Eq.  (1),  then  we 

obtain: 

"ae  =  ?(T^e>T^He^0e>-{0e>[Ge]{5*}  +  {Qo}T{5e})  (5) 

0  0 

where  [K  ] ,  [G  ]  and  {Q Q }  are  combined,  respectively,  according  to 
the  specified  sites  by: 
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CH-1  J  =  JVni[A‘i’  ]T[S<i>  ][A'i’]dv 

CG,i,D  =  Jdvo.[R-.-:TCL‘i>  Ids 

{Q  i,}T  =  -fSa  .<  Pli’  }T[L  »:ds 

ouni 

0  0 

(3  }  and  {6  }  are  stress  parameters  and  nodal  point  displacement 
vector  (in  the  broad  sense)  of  the  entire  element,  and  combined, 
respectively,  with  {3^}  and  {6^}  (i=l,  2,  ...,  m)  . 


Stresses  o  ,  a  and  a  should  continue  among  layers;  this 
xz  y  z  z 

continuation  condition  can  be  expressed  as: 


CctxV’  ,  a y?  ,  o2“’  ]  =  [  0  ,  0  ,  0  3 
.  °ri* "  -  °z(U  M  3-  [Ox'i*  ,  a#  ,  a2“’  3  =  [  0  ,  0  ,  0  3 
Cox'f’  -  cr7‘?>  ,  a2<mi  3  =  C  0 , 0  ,  0  3 


We  substitute  Eq .  (2)  into  (6) ;  this  continuation  condition 

can  be  expressed  with  the  stress  parameter  as: 

LBjm  =  {0>  (7) 

Using  Lagrange's  operator  on  Eq .  (7)  and  introducing  into 

arbitrary  function  (5) ,  then  the  following  equations  are  obtained 
through  differentiation  and  rearrangement: 

nme=  E(l{6e}T  [K03{5*}  -  {Q0}T{6e})  (8) 

{Pe}  =  (CHe3'1[Ge3-[He3',[B3T([B3[H'3'1CB3T)*1[B3CHe3',[Ge3){5e}  (9) 

In  the  equations, 

CK03  =  [Ge3[He3*1[G'3-[Ge3[He3"![B3T([B3[He3'1[B3T’l[B3[He3'1[Ge3 

This  is  apparent  from  the  above-mentioned  deduction  process. 
Since  stress  {cr^}  is  different  for  each  layer,  and  the  elasti¬ 
city  matrix  [s^]  of  each  layer  is  generally  different,  therefore 
the  stress  field  {e  }  =[s  ]{a  }  of  each  layer  differs  layer  by 

layer.  Based  on  geometric  equations,  it  is  known  that  there  are 
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different  variation  (in  terms  of  coordinates)  equations  for  the  dis¬ 
placement  of  each  layer  with  the  exception  of  transverse  displace¬ 
ment  w;  this  is  the  different  function  of  the  coordinates.  Hence, 
the  post-deformation  normal  line  is  partitioned  on  different  seg¬ 
ments  according  to  the  layer.  This  element  can  be  used  to  take  into 
account  the  local  warping  effect. 

On  the  element  boundary,  according  to  the  hypothesis  of  YNS 
theory,  we  have: 

{&'}=Cw]{qf}  (10) 

where  {q  }  is  the  nodal  point  displacement  (in  the  broad  sense)  at 
the  midplane  of  the  laminated  plate. 

Substituting  Eq .  (10)  into  (8),  we  obtain: 

*me=  {  qe}T[Ke]{qe}-  {Qe}T{q*}  (11) 

In  the  equation,  [Ke]  =  [W]T[K0 ]'W]  is  a  single-element  rigidity 
matrix:  {Qe}  =  {Q0 }T[W]  is  the  single-element  nodal  point  as  vector. 

From  dynamic  analysis,  the  modified  variation  principle  is: 

=  I"  ?  {|,[|  7.i  1 W"  >Te  W“  Wv.  |  Vn.  i<a  'i'  >T[s  <f  ](0  <P  >dv- 

r  -  (12) 

-Javn,  <p '!>T{u  + 

where  is  the  mass  density  of  the  i-th  layer;  { u ^  }  is  the 

derivative  with  respect  to  time  of  the  internal  displacement  vec¬ 
tor  of  the  i-th  layer. 


In  Eq .  (12),  in  principle  the  internal  displacement  field 

{uj}  and  single-element  boundary  displacement  field  can  be  inde¬ 
pendently  assumed.  However,  for  convenience  the  paper  adopts  the 


internal  displacement  field  compatible  with  the  boundary  displace¬ 
ment  field.  We  assume  that  the  internal  displacement  field  is: 

{uj"1  }  =  [X  }  (13) 

Substituting  Eqs .  (2),  (3),  (4)  and  (13)  into  Eq .  (12),  we 

obtain  the  following  after  rearrangement: 

f”c(l(q(TCM-]<q}-4<,*}CK']{q'})dt  (14) 

J  t1  a  1  Z 

0 

In  the  equation,  [M  ]  is  the  mass  matrix  of  a  single  element  and  is 
composed  of  the  mass  matrices  of  all  the  layers.  When  the  free 
vibration  problem  is  considered,  no  load  item  appears  in  Eq .  (14). 

Jvn.p'i>CN';>]T[N(i>  jdv 

From  Eq.  (14)  a  characteristic  equation  with  the  intrinsic 
property  of  power  can  be  derived: 

(CKJ-olCM3){q  }  =  Q  (15) 

In  the  equation,  [K]  is  the  overall  rigidity  matrix;  [M]  is  the 
overall  mass  matrix. 

III.  Hypotheses  of  Field  Variates 

(1)  In  a  single  element,  the  stress  field  hypothesis  of 
each  layer  is: 

<*x>  =  P  +  04u'  X  +  P  «>  y  +  3  z  +  3  xz  +  P  y  z 

=  3*"’  +  P,"’  x  +  3,“’  y  +  ^  <V  z  +  3  ,‘V  xz  +  3,‘V  yz 

=  3,“’  +  3«(i'  x  +  3.“’  y  +  0  ,'**  z  +  3,‘V  xz  +  P  <V  yz 

2  (16) 

° j*  =  -  3 *  -  p  *  -  y  o  ,‘i*  +  p  ;v ) + p  ,‘i’ 

=  -  p.*”  x  -  p.‘"  z  -  y  (P  iv  +  a  /i’ ) + p  ,‘v 

a*'”  =  0 
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In  the  equation,  z  is  the  distance  from  a  point  in  each  layer  to 
the  midplane  of  that  layer.  Equation  (16)  satisfies  the  stress 
homogeneous  equilibrium  equation. 

(2)  The  boundary  displacement  field  of  each  layer  in  a 
single  element: 


We  take  displacements  u,  v  and  w  along  the  interlayer  direc¬ 
tions  x,  y  and  z,  considered  as  the  nodal  point  displacements  (as 
in  Fig.  1)  in  the  broad  sense.  Here,  the  author  and  his  colleagues 
assume  that  w  is  a  constant  along  the  entire  thickness  direction 
of  the  single  element  in  order  to  be  consistent  with  c  =0.  For 
boundary  AB,  the  equation  expressing  boundary  displacement  is: 


“A'i  =  j(u  +U5''‘  )(l  - —)  +  i-Cu,11’  + 


+  iTirC(u.<i‘  + 


-  =  +  v  j’’  )(1  -  +  y(v;  ‘  1  +v,‘i>)-f-+h~ir[(v.'i 


a 

x  z 


(17: 


-u;ii)d--y)+  (v,p  - vj*’ 


Wab‘  =  »l"'  (I  -  ~)  +  w,"1  -^-=  W,(l  -  -J-)  +  ~ 


'  1> 


In  the  equation,  h^  is  the  thickness  of  the  i-th  layer.  Displace¬ 
ments  of  other  boundaries  can  be  written  out  following  the  previous 
line  of  reasoning. 


(3)  Transformation  matrix  [W] : 

We  take  four  corner  points  at  the  geometric  midplane  (of  the 
laminated  plate)  as  nodal  points;  there  are  five  degrees  of  free¬ 
dom,  three  displacement  degrees  of  freedom  (u,  v  and  w) ,  and  two 
degrees  of  freedom  of  the  angle  of  rotation  (namely  the  angle  of 

rotation  6  around  the  x-axis  and  the  angle  of  rotation  9  around 
x  y 

the  y-axis)  for  each  nodal  point,  as: 
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(18) 


“i'*’  =  U,  +  ljBy 
v  -  V)  r  1;9X| 
w,‘  =  w. 

In  the  equation,  1^  is  the  distance  from  the  point  (to  be  solved) 
of  nodal  point  displacement  to  the  geometric  midplane.  The  other 
nodal  point  displacements  of  each  layer  in  the  broad  sense  can  be 
similarly  expressed  in  terms  of  the  nodal  point  displacement  on 
the  geometric  midplane.  From  Eq .  (18),  the  transformation  matrix 

[W]  can  be  formed. 


Fig.  1.  Nodal  point  displacement. 


(4)  Hypothesis  of  internal  displacement  field: 


*y 


y,li)  =d- 


f  +  Tb->*  +  o-r-f +  -ib-K*  +  h08- 


(i9: 


+  (- 


xl  )v,+(  x  -  ,  xy 

ab  a  a 


b  )(2+hi)0*=  +  A/b  v3+~y-(Z  +  hi)0XJ  + 
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+  <-6  - 


iv-)  v  +(>'. 

ab  '  4  ^  b 


xy 

ab 


)(z  -  h,)Ox* 


(19) 


(1- 


X 

a 


xy 

T_ab')Wl 


-L-)w 
a  b 


ab 


(V 


x  y_  x 
ab  ^  w  * 


where  hi  is  the  distance  from  the  midplane  of  the  i-th  layer  to 
the  geometric  midplane  of  the  plate.  It  can  be  proved  that  Eq . 
(19)  is  compatible  with  the  displacement  field  (of  each  boundary) 
at  the  boundary. 

IV.  Numerical  Results 

(1)  Flexural  stress  analysis  of  simply  supported  rectangular 
laminated  plate  [0°/90°/0°] 

The  plate  dimensions  are  3axa;  the  bearing  load  is 

q  =  q.Sir  jr^in~3a~  .  The  origin  of  the  corresponding  coordinate 
system  is  taken  at  a  corner  point  of  the  plate;  the  y-axis  extends 
along  the  long  side;  the  x-axis  along  the  short  side;  and  the 
ratio  between  plate  width  and  plate  thickness  is  (a/h)=10.  In 
[3],  Pagano  gives  an  exact  solution  for  this  problem.  In  [1] 
a  finite  element  analysis  is  made  for  this  problem;  the  numerical 
results  and  the  calculated  results  in  the  paper  are  listed  in 
Table  1. 


Table  1.  Stress  and  deflection  of  laminated  plate. 


( 

a)  * 

•i 

(bP 

IS- 

Tc T~ 

a 

3< 

a. 

Ox. 

O  y  z 

CJ  x  y  W 

dtj 

!  A 

<-v- 

3a 

2 

2  V2 

3a 

2 

4>». 

-T"0,(f 

,o.o)(o,o.  ± 

(d) 

Cl  )  4M. 

2  X  1 

54 

±  0 

828 

±  0.0462 

0.287 

0.0104 

±0.0136  1 

0131 

(e. 

JM. 

2  X  1 

30 

±  0 

834 

l 0.0467 

0.283 

0.0101 

±0.0139  1 

0133 

(d) 

(!)♦<, 

4x2 

135 

±  0 

74  1 

a  0.0433 

0.386 

0.0139 

±0.0126  0 

9427 

(e) 

*  *.*  *. 

4x2 

75 

t  0 

744 

± 0.0440 

0.383 

0,0136 

±0.0128  0 

9430 

(d) 

Cl)**. 

8x4 

405 

±  0 

714 

t  0.0426 

0.415 

0.0150 

±0.0122  0 

9255 

(o) 

8x4 

225 

±  0 

737 

±  0.0428 

0.413 

0.0149 

±0.0123  0 

9260 

(f ; 

in  m 

0 

-  0 

726 

725 

0.0413 
-  0.0435 

0.420 

0.0152 

-0.0120  Q 
0.0123 

919 

(Key  on  following  page.) 
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Key  to  Table  1:  (a)  Method;  (b)  network;  (c)  degree  of  freedom; 

(d)  single  element  from  [1]?  (e)  single  element 

in  the  paper;  (f)  exact  solution 

Table  2.  Effect  (^=10)  of  dimensionless  frequency 
2  =  oa:(p/E;h:)~-~  of  square  plate  by  layer  number. 


a)  &  it 

(b)C 

(c)  a  C  5  ) 

(d)  $.  X. 

2 

1 5.714 

13.04 

13.581 

4 

13.060 

18.46 

18.385 

6 

18.295 

19.09 

19.003 

8 

19.023 

19.23 

19.230 

10 

to 

o 

- 1 

4~ 

19.38 

19.351 

1  2 

10.003 

— 

19.429 

14 

13.113 

— 

19.484 

18 

19.122 

18.48 

19.500 

Key:  (a)  Layer  number;  (b)  single  element  from  [4];  (c)  ref¬ 
erence  [5];  (d)  this  paper. 


(2)  Intrinsic  property  analysis  of  bias- laid  laminated  plate 
The  fiber-laying  angle  of  the  plate  is  [-4 5 °/4 5 °/-45 °/4 5 °/ . . . ] ; 
the  elastic  constant  of  the  material  is 

G  1  3  G;  3 


Ei  .  G ! z 
E  =  40,  £ ,  -  0  <  6, 


=  0.5,  v  =  0.25. 


The  boundary  condition  of  the  plate  is  u=w=0,  0=0.  At  a, 

A 

where  x=0,  v=w=0,  and  0^=0.  At  b,  the  dimensionless  frequency 
values  of  a  plate  with  different  layer  numbers  are  listed  in 
Table  2. 


Va  Conclusion 

From  the  above  analyses  and  numerical  calculations,  the 
following  conclusion  can  be  obtained:  to  conform  to  the  properties 
of  a  laminated  plate,  stress  parameters  should  be  assumed  layer  by 
layer.  However,  at  the  boundary  of  a  single  element,  the  displace¬ 
ment  hypothesis  in  the  YNS  theory  can  be  adopted.  A  single 
element  thus  constituted  can  not  only  take  into  account  the  effect 
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of  shear  deformation  and  the  local  warping  effect,  but  also 
involves  fewer  degrees  of  freedom  of  a  single  element. 
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